{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "from DFT_WW import *"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "nteract": {
     "transient": {
      "deleting": false
     }
    }
   },
   "source": [
    "The DOS and XPS data can be imported however you choose. Here, we use np.loadtxt to import from CSV files. By default, this requires the file to have no headers. The density of states can be unphysically jagged; if so, we recommend applying some smoothing."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false,
     "source_hidden": false
    },
    "nteract": {
     "transient": {
      "deleting": false
     }
    },
    "outputExpanded": false
   },
   "outputs": [],
   "source": [
    "dosFile = 'Pt_DOS.csv'\n",
    "xpsFile = 'Pt_XPS.csv'\n",
    "\n",
    "energyGrid, dos = np.loadtxt(dosFile,delimiter=\",\",unpack=True)\n",
    "energyXPS, XPS = np.loadtxt(xpsFile,delimiter=\",\",unpack=True)\n",
    "paramsDict = {'bindingEnergy':0.26,'intensity':1,'lorentzianWidth':0.48,'gaussianWidth':0.62,'scaleDOS':0.00050625}\n",
    "name = 'DFT+WW'\n",
    "fermiLevel = 0\n",
    "\n",
    "visualizeInput = False\n",
    "visualizeOutputs = False"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false,
     "source_hidden": false
    },
    "nteract": {
     "transient": {
      "deleting": false
     }
    },
    "outputExpanded": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEaCAYAAAAVJPDdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3dd3yV9dnH8c+VHZKQEAgEEkbYW4UAERc4cRWr1rpHXTiqtmr1aWu1+rTWttqn1lVcdeNWVBQVFZEhG2QTwgojCSshITvX88e5sSchITkk59wnOdf79cor55x7fc9NONe5x+/3E1XFGGOMOSjM7QDGGGOCixUGY4wxtVhhMMYYU4sVBmOMMbVYYTDGGFOLFQZjjDG1WGEwbY6InCAia93O0RQisklETnU7hzHerDCYFuN8yJWKSLHXzxOBzqGqs1R1QKC3KyK9RERFJKKF1veuiEyu89oHB/epiFwtItXOfi4SkaUico7XvL8VkY3O9FwRefMw2zr4b7dfRPaJyBwRmSQiTfqMaOn3btxlhcG0tHNVNd7r59ZAbryNfTDdAlwgIuMBROTnwDHAvV7zzFXVeCAJeB54S0SSReQq4ArgVGd6JjCjke2dq6oJQE/gL8A9zjpNiLHCYAJCRJ4WkXe8nj8iIjPEY5zzjfa3IrLL+fZ6mde80SLydxHZIiJ5IvKMiMQ60w4ue4+I7ARePPia1/KbRORuEVkuIiUi8ryIdBGRT51vyF+KSAev+bOcb8z7RGSZiIzzmvaNiDwkIrOdZT8XkU7O5G+d3/ucb+nHikgfEflKRHY77+01EUlqyj5T1Z3AncCzItIDeBy4UVWL65m3BngBiAV6A6OA6aq64eC6VHVy3eUa2G6hqk4Ffg5cJSJDnfd+togscY5OtorIA16Lteh7N+6ywmAC5U5guHP64wTgWuAq/W+fLKlAJyANuAqYLCIHTwc9AvQHjgb6OvP8wWvdqUAynm+6NzSw/QuA05z1nAt8CvzW2WYYcBuAiKQBnwD/66zzLuBdEUnxWtelwDVAZyDKmQfgROd3knO0NBcQ4GGgGzAI6A48cPhd9V+q+h9gA7AY+ExVP6tvPudI6TqgGFgPzAOudApipoiEN3WbXtueD+QCJzgvlQBX4jk6ORu4SUTOc6a1+Hs37rHCYFraB8437YM/1wOo6gHgcuAx4FXgl6qaW2fZ+1S1XFVn4vlwvkhEBLge+JWq7lHV/cCfgYu9lqsB7neWLW0g179UNU9VtwGzgO9VdYmqlgPv4zlFg5NxmqpOU9UaVf0CWAic5bWuF1V1nbOtt/AUrHqparaqfuFkK3De/0mH2X/1mQV0xLPf6soSkX3ATuAS4KfON/5XgV8CZwAzgXwRubee5RuzHU+BRFW/UdUfnP2yHHjjcO+lhd67cUFbOh9rgsN5qvplfRNUdb6I5OD5pv1Wncl7VbXE6/lmPN80U4B2wCJPjQA830S9vwEXqGpZI7nyvB6X1vM83nncE/iZiJzrNT0S+Nrr+U6vxwe8lj2EiHTGcwroBCABz5exvY1k9V6+H54jkqeAR0VkpKpWes0yT1WPr29ZVX0NeE1EIoHznMdLVHV6U7eP5+hsj5NlDJ5rD0PxHClFA28fJnuz3rtxjx0xmIARkVvwfJhsB35TZ3IHEYnzet7DmW8Xng/uIaqa5PwkOhdUD2rJLoK3Aq94bStJVeNU9S9NWLa+HA87rw9X1fZ4jkiknvkO4RwtPQf8H55v/yV4Lgj7RFUrVfVtYDmeD/UmEZFReArDd85LrwNTge6qmgg8w3/fS4u+d+MuKwwmIESkP57z9pfjuVvmNyJS9xTMH0UkyrkGcQ7wtnNR9VngH843UEQkTUTO8FPUV4FzReQMEQkXkRjnYnZ6E5YtwHNaq7fXawl4zvvvc65f3O1DlpvwXAP5s7MfrsWz3wY2tqBzLedsEUkQkTARORMYAnzfhGXbi+e21ynAq6r6g9d72aOqZSIyGs+1loNa+r0bF1lhMC3tI6ndjuF958Loq8AjqrpMVdfjufD7iohEO8vtxHOaYTvwGjBJVdc40+4BsoF5IlIEfAn4pZ2Cqm4FJjr5CvAcQdxNE/6vONdR/gTMdq6vZAF/BEYAhXium7zXlBwi0h3PtZRrVbXCWf8q4FE8dyk19s27yHkPW4B9wF+Bm1T1u8Ms85GI7Mfznn+H55rANV7TbwYedOb5A16nA1vyvRv3iQ3UY9zm3A76qqo25Vu5McbP7IjBGGNMLVYYjDHG1GKnkowxxtRiRwzGGGNqafUN3Dp16qS9evVyO4YxxrQqixYt2qWqKfVNa/WFoVevXixcuNDtGMYY06qIyOaGptmpJGOMMbVYYTDGGFOLFQZjjDG1WGEwxhhTixUGY4wxtQSsMIjICyKSLyIrGpguIvK4iGQ7QzCOCFQ2Y4wx/xXII4b/ABMOM/1MoJ/zcwPwdAAyGWOMqSNghUFVv8UZCaoBE4GX1WMekCQiXQOTzphDHaiooqyy+sfny7buY8W2Qsoqq1m1vYjXv9/CruJyFxMa4x/B1MAtDU8/8AflOq/tqDujiNyAM+h7jx49AhLOhJ7LnvuemhrlzRuPZVdxOROfnH3IPP+csY7HLjqa4/p2ciGhMf4RTIWhvoFH6u3hT1UnA5MBMjMzrRdA0+J2FZezZMs+AH7/wQp2FnqGlJ54dDfSkmJR4NjeHbnz7WVc9tz3PHtlJqcN7uJiYmNaTjAVhlygu9fzdDyjeRkTcHM37Abg1EFdeGdRLgB//MkQrhrbq9Z83949nhP++hVTl223wmDajGAqDFOBW0VkCjAGKFTVQ04jGeNPZZXVfLBkG8/OyiEhOoKnLx/B9JU72V1cwRVZPQ+ZPzYqnFMHdeHj5TvIKyojOS6K9xbnMmFoVxJjI114B8Y0X8AKg4i8AYwDOolILnA/EAmgqs8A04Cz8Izte4DaY80a43eV1TVc9cJ8vt/ouUfiofOGEhkexjnDux12uYtGdee9JdsY8+cZP762YlsRD5031K95jfGXVj9QT2ZmplrvqqYlPPHVev7++TruPmMAR3dP8umC8sZdJUxZsIV/z8whKiIMVc9F6xE9OvgxsTFHTkQWqWpmvdOsMJhQV1FVw5sLtvDAR6s4c2gqT1x65G0rcwqKiYoI47LnvqewtJLPbj+R1MSYFkxrTMs4XGGwLjFMyFJVvl6Tz8+emcN9H65kSLf2/OWC4c1aZ++UeNI7tOPFq0dRXlnD+U/NZl7O7hZKbExgWGEwIUlV+cOHK7nmPwtYllvIr0/rz4e3HEd8dMtcduudEs/zV2VSrcrjM9a3yDqNCZRguivJmIBQVf740SpembeZK7J68vNR3RnSrT0i9TWlOXJj+3biwpHpPDMzh70lFXSIi2rR9RvjL3bEYELO899t5D9zNnH9CRk8OHEIQ9MSW7woHHTGkFSqa5QvV+f5Zf3G+IMVBhNStu0r5dHP13HqoM789qxBfisIBw1LS6RbYgzTV1phMK2HFQYTUv755TqqVXngJ0P8XhQARITTh6Qya30BJeVVft+eMS3BCoMJGXlFZbyzKJfLx/QkvUO7gG333KO6UV5Vw1sLtzY+szFBwAqDCRkfLdtOjcJlWYHtkXdkzw6Mzkjm3zNz7KjBtApWGExIWLW9iBdnb2JoWnv6pMQHfPv3TBjAzqIyhtw/nSe/zg749o3xhRUG0+ZV1yi3vrGY8qoa7jt7sCsZRvZM5q7T+wPwnzmbqKqucSWHMU1hhcG0eW8v3EpOQQkPThzCmN4dXctx68n9ePqyERTsL2fOBmsNbYKXNXAzbdbD01YTHia8MHsjYzKSmTAk1e1IjB/YmaiIML5dV8CJ/VPcjmNMvawwmDZpwaY9/PvbHABS28fwxKUjCAvz/+2pjYmJDGdEjyTmWv9JJojZqSTTJr00ZxNxUeFckdWT56/OJCUh2u1IPzq2dydW7ShiT0mF21GMqZcVBtMmLcvdx7iBnXnovKEM6ZbodpxaJgxNRcA61zNBywqDaXMKSyvZuqeUId3aux2lXgNSE7g8qycvz93Eim2Fbscx5hBWGEybs2p7EUDQHSl4u/P0ASTHRfGnT1a7HcWYQ1hhMG3OfGfM5sFdg/OIASAxNpJfHJ/B3Jzd5BQUux3HmFqsMJg2Jb+ojMnfbmD8gJSguuBcnwtGpBMm8M6iXLejGFOLFQbTpkxdtp2Simp+f447LZx90aV9DOMGdObdxbnWEtoEFSsMps14fMZ6/veT1fTrHO9Kf0hH4qLMdPKKypm1fpfbUYz5kRUG0+qt2FbIHVOW8NgX6wA4oV/raVF88sAuJMdF8fYi65LbBA9r+WxatfKqaq56YT67ncZit53Sj6vH9nI3lA+iIsL46TFpvDx3E3tKKki2caFNELAjBtOqfbxsB7tLKsjqncwdp/bj16f1b3Ufrj89Jo3KamWGjQttgoQdMZhW7YOl2+ieHMsb12cFZKhOfxjSrT2p7WP4ak0+P8vs7nYcY+yIwbRe+fvLmJ29i4lHpbXaogCecaFPHtSZb9YWkJ2/3+04xlhhMK3XR8t2UKNw3jHd3I7SbDeP60NcdAS/fW+F21GMscJgWqeaGuXdRbkM6daevp0T3I7TbOkd2nH2sFRW7ShCVd2OY0KcFQbTKr0ybzOrdhRxzXEZbkdpMb1T4ikur6Jgf7nbUUyIs8JgWp2yymr+9dV6juvbkQtGpLkdp8X0TokDYENBictJTKgLaGEQkQkislZEskXk3nqmJ4rIRyKyTERWisg1gcxnWoepS7ezq7iCm8f1bdUXnevq7bTWztllneoZdwWsMIhIOPAkcCYwGLhEROp2aHMLsEpVjwLGAY+KSOu6Kd34lary3Hc5DExNYGyfjm7HaVFd28cQGxnOrHW7qKmx6wzGPYE8YhgNZKtqjqpWAFOAiXXmUSBBPF8D44E9QFUAM5ogVlJexXlPzWFdXjHXndC7TR0tAISFCdcen8FnK3fym3eXW3EwrglkYUgDvDuEyXVe8/YEMAjYDvwA3K6qh3Q7KSI3iMhCEVlYUFDgr7wmyHy1Jp9lW/dxRVZPzju69d+iWp87T+/PL0/uyzuLcvnSWkIblwSyMNT39a7uV6IzgKVAN+Bo4AkROWS0FVWdrKqZqpqZktJ6OkwzzfPZyp10io/mgZ8MISK8bd43ISLcdko/OsVH8+5iG6fBuCOQ/7tyAe/2/ul4jgy8XQO8px7ZwEZgYIDymSCWnb+fL1blcebQVMLD2tYppLoiw8M47+hufLUmn/1llW7HMSEokIVhAdBPRDKcC8oXA1PrzLMFOAVARLoAA4CcAGY0Qai0oprbpywlLiqc207p53acgDh5YGcqq/XHYUqNCaSAFQZVrQJuBaYDq4G3VHWliEwSkUnObA8BY0XkB2AGcI+q2ggmIe6Rz9awakcRj150VNAP19lSRvTsQExkGN9l25+/CbyA9q6qqtOAaXVee8br8Xbg9EBmMsEtp6CYV+Zt5pLRPTh5YBe34wRMTGQ4o3olMyd7t9tRTAhqm1fwTJvx6rwthAn86tT+bkcJuMyeyazL30+RXWcwAWaFwQSt8qpqPli6jdMGdwmZU0jejumRhCos31rodhQTYqwwmKD19sJc9pRUcOnonm5HccXRPZIQgUWb97odxYQYKwwmKO0uLudfX61nRI8kjuvbtrq+aKr2MZGM6NGBKQu2UFpR7XYcE0KsMJig9Nv3f2DvgUoenDi0zXV94Yt7JgxkR2EZ71hjNxNAVhhM0Fmzs4jpK/O4ZVxfhqYluh3HVaMzkumeHMt3663rFxM4VhhM0Hnmmw3ERYVz9dhebkcJCmN7d2Jezh6qrVM9EyBWGExQ2brnAB8t38GlY3qQ2C7S7ThBYWzfjhSWVrIsd5/bUUyIsMJggsqzs3IIE7j2+N5uRwka4wd2JiYyjLcXbm18ZmNagBUGExRqapSX527izQVbOf+YdFITY9yOFDTax0RyzvBuTF26neJyG57E+J8VBhMUfvPucv7w4UqOSk/iV6eFXivnxlwyugclFdV8tKxuh8TGtDwrDMZ1Czft4Z1Fudx4Ym/evDHLjhbqMaJHEgO6JPDG/C1uRzEhwAqDcd3r87eQ1C6S20/tF9JtFg5HRLhkdHeW5xayYpt1kWH8ywqDcd38jXvIyuhIu6iAdvbb6vz0mHSiI8J4yy5CGz+zwmBctaOwlNy9pYzKSHY7StBLbBfJSf1TmLE6H1Vr02D8xwqDcU3u3gNM+2EnAKN7WWFoivEDO7NtXynr8ordjmLasEaP3UWkKf9ja1TVWt+YJttfVsnxj3wNwNC09gxNa+9yotZh/IDOAHy5Oo8BqQkupzFtVVNO6m53fg53VTAc6NEiiUxI+HJ13o+P7zt7sF10bqLUxBhG9Eji4+U7uGV8X7fjmDaqKYVhtaoec7gZRGRJC+UxIeLDpdtJS4pl1m/GExZmRcEX5wzvxoMfr2JDQTF9UuLdjmPaoKZcYzi2heYxhsLSSh7+dDXfrC3g/BFpVhSOwFnDuiICHy/b4XYU00Y1WhhUtawl5jEG4O/T1/LvmTkAXDYmNEdma67UxBhG9Uzm4+XWCtr4h92VZAImr6iMNxdu5aj0RF69doy1cG6Gs4alsj6/mI27StyOYtqgZhUGEUltqSCm7Xv22xyqa5R/XTKC4/t1cjtOqzZ+oOfupJlr811OYtqi5h4xPN8iKUyb933Obl6et5mJR3ejR8d2bsdp9Xp2jCOjUxwz19nIbqblNaswqOrZLRXEtF2rthdx3UsL6d4hlt+fPdjtOG3GSf1TmJuzm7LKarejmDbGrjEYv/vztNXERIXz6nVjSI6LcjtOm3HSgBTKKmuYv3GP21FMG9PkwiAiC0TkeRG5Q0ROFpEUfwYzbcPm3SV8l72LK7N60jUx1u04bUpWRkeiIsL4YlVe4zMb4wNfjhgmAm8DUcAkYJOIbPZLKtNmPDsrh8hw4WeZ3d2O0ubERoXzk6O6MWXBFnIKrO8k03KaXBhUdbuqfqaqf1XVi4BM4Dn/RTOt3dY9B3hzwVZ+Pqq73ZrqJ/dMGIiI8Nr3NoCPaTm+nEqq1ReSqq4GhrR4ItNmPD5jPSLCreP7uR2lzUpJiGZMRjLf2G2rpgX5cirpTRHJFZFZIvKUiDwGDPRlYyIyQUTWiki2iNzbwDzjRGSpiKwUkZm+rN8Ejz0lFby3ZBuXju5hRwt+Nm5AZzYUlLB1zwG3o5g2oslDZqnqsQAi0hcYBiQDjzV1eREJB54ETgNygQUiMlVVV3nNkwQ8BUxQ1S0i0rmp6zfB5cvVeVTXKBeMSHc7Sps3tk9HABZt3kv3ZGsjYprP57EUVTUbyD6CbY0GslU1B0BEpuC5oL3Ka55LgfdUdYuzLTs+boUe+2Idj89YT7fEGBtnIQD6dY4nOiKMldsLOe+YNLfjmDYgkO0Y0gDvwWpznde89Qc6iMg3IrJIRK6sb0UicoOILBSRhQUF1vIzmMzL2f1jUbj3rEE2zkIARISHMTA1gRXbityOYtqIZo2+LiKpqrqzqbPX81rdgWsjgJHAKUAsMFdE5qnquloLqU4GJgNkZmba4LdBQlX562dr6JYYw4w7xxEbFe52pJAxJC2Rj5Ztp6ZGrStz02yB7CspF/C+mT0dz8hwdef5TFVLVHUX8C1wVPMimkCZ9sNOFm/Zxy0n97WiEGDH9enE/rIq3ly4tfGZjWlEIPtKWgD0E5EMEYkCLgam1pnnQ+AEEYkQkXbAGGB1czKawPhmbT6/ff8HhqUl8nNrzBZwZw1LZXRGMo9+vo6q6hq345hWzpd2DB1F5CYRuUZERouIT/0bqGoVcCswHc+H/VuqulJEJonIJGee1cBnwHJgPvCcqq7wZTsm8GpqlLvfWU5KQjRPXHoMEeHWBVegiQjXjO3FruJyvre+k0wz+XKN4X3gS+AmYB1wrIjkqGqT2zKo6jRgWp3Xnqnz/G/A33zIZVy2akcRBfvLuWfCQHp2jHM7TsgaP7AzcVHhfLRsO8f1tfEuzJHz5atdgqo+COSp6knAJcCL/ollWpMvV3s6cTupv/Wr6KaYyHBOH5LKpyt2UlFlp5PMkfOlMBwc17lcRGJV9V3gLD9kMq3Iurz9/HtmDif2TyElIdrtOCHvnOFdKSyt5Ltsu43bHDlfTiX9XUSSgTeBF0RkDoe2QzAhYsbqPJ7+ZgOLt+ylY3w0f7twuNuRDHBCvxTax0Tw0bIdnDywi9txTCvlS++q76rqHlV9DM91gu54Wi6bEFNaUc3tU5ZSUFzOKYO68Mq1o+nS3vpDCgZREWGcObQrX6zKs5HdzBE7ogZuqvpKSwcxrcfnq3ZSXF7F5CtHMraPXeQMNuce1Y03F27lm7X5TBja1e04phVq9IhBRBa3xDym7Xhl7mbSkmLJyujodhRTj6zeyXSKj+KjZTvcjmJaqaYcMQwSkeWHmS5AYgvlMUGsrLKaP32ymoWb9/LAuYOt64UgFREexlnDuvLWwq0Ul1cRH92snm9MCGrKX0xT2inYycwQ8OLsTbwybzNjMpL5+agejS9gXHPuUd14ee5mZqzOY+LRdo+I8U2jhUFVbVxnw87CMiZ/u4ET+6fw8i9Gux3HNGJkjw6kto/ho2XbrTAYn1nfBaZRqspNry2islr53VmD3I5jmiAsTDhneFdmriug8ECl23FMK2OFwTRqxup8lmzZx33nDGJAaoLbcUwTnXtUNyqrlemrmtozvjEevnSi96WIWBfYIaa8qpqHP11Nj+R2nG/DdLYqw9MT6Z4cy+cr89yOYloZX44YfgP8Q0ReFBG7OToErN5RxBXPzWdDQQkPThxCpPWa2qqICKN6JbN06z5UbTwr03S+tHxerKonAx8Dn4nI/b52vW1aj+z8YiY+OZv1+fv564XDGTegs9uRzBE4Kj2JXcXl7Cgsa3xmYxw+fQUUzwC+a4GngV8C60XkCn8EM+4pr6rmf95bTkxEGJ/dcSIX2cA7rdbwdE8To+W5+1xOYloTX64xfAdsA/6Bp/O8q4FxwGgRmeyPcCbw9pRUMOmVRSzYtJcHJw61PpBauUFd29MuKtyuMxif+NIkchKwUg89WflLEbHhN9uAFdsKufCZOVRWK3/+6TDOO8buf2/tYiLDuXhUD16eu4k7zxhAWpKd/TWN8+VU0st1i4KIZDkPfRn72QSpF2ZvJCIsjGm3ncClY6xlc1tx3QkZADw3K8flJKa1aEoneheJyF+ABBEZJCLhXpMnA6iq/cW1coWllXy8fAfnHdPN2iq0Md2SYvnJ0d2YMn8re0sq3I5jWoGmHDHMBlYBHYDH8FxwXiwiHwOl/gxnAufzlZ7hIC8caRea26JJJ/WhtLKal+ZucjuKaQWa0lfSNuBlEdmgqrMBnJHcMoA1fs5nAiCnoJjHv1pPeodYjkq3jnLbov5dEjh5YGemzN/KHaf2dzuOCXJNOZV0v/NQRSQBwBnJbZGqlvg1nfG70opqrnxhPvlF5dw6vi+eO5JNW3R8307sLCojf7+1aTCH15S7kj53ft8BDBWRCDynlpYDy1X1HX+FM/737KwccveWMuWGLLJ628A7bdmQbu0BWLm9iM4D7DZk07BGjxhUda7z+yJVHQwMAx4ENgBZh1vWBLfs/GJemL2RUwd1tqIQAgYfLAzbCl1OYoKdLw3cZopIe1UtB0bjuRj9W78lM3513wcrOPWxmRwor+bWk/u5HccEQEJMJBmd4vh2/S7rO8kcli/tGJJUtUhERgLX4ykMz/onlvGXmhplXd5+Xpm3mQtHpjPzN+M4unuS27FMgFx5bE/mb9zD9JXWFbdpmC+FodK5vnAl8Iiq3g8M8U8s4w8rthVy+v99y+n/+BaAu88YQNdEawkbSq7I6snA1AQe+ng1pRU2Iq+pny+F4V/AMuAc4CPntfgWT2T8Yk9JBZc8O4+S8irOPaob15+QYf0ghaCI8DDuO2cw2/aV8rkN4GMa0OS+klT1JRF5F6hW1VIR6QvM9V8001KqqmuYunQb+8uqeP26LIZZW4WQNiYjmXZR4SzZss/Ggzb1anJhEJGTgcuAfSKyAs/tqpP8Fcy0jNy9B/jFfxawLq+YHsntGJrW3u1IxmUR4WEMT09kyZa9bkcxQcqXU0mv4hmkZx7QG/gDsNIfoUzLeXjaGrbtLeWk/incfko/a8BmABjRowMrtxdRUl7ldhQThHwpDNmq+r6qvq2q96nqRFXt68vGRGSCiKwVkWwRufcw840SkWoRudCX9Zv/Wp+3n+Mf+YpPftjB5Vk9eekXo7lgpI3ZbDxOGdSFqhrl9e+3uB3FBCFfCsNMEfmVHOFXTqdX1ieBM4HBwCUiMriB+R4Bph/JdgxUVNVw9YsLyN1bSp+UOK4a28vtSCbIjOzZgeP6dmTyrByqqmvcjmOCjC+FYQhwE7BDRD4RkT+JyM98WH40nqOOHFWtAKYAE+uZ75fAu0C+D+s2Xj5dsYNt+0p58epRzLhzHN1scBZTjyuP7UXB/nJmZe9yO4oJMk0uDKp6vqr2x9Or6v3AenzrEiMN2Or1PNd57Ucikgb8FHjmcCsSkRtEZKGILCwoKPAhQtu3p6SCv01fS5+UOE7qn+J2HBPExg/oTFK7SN6cv7XxmU1I8eWupH7AvUCpqt4KLPRxW/WdgqrbLv//gHtUtfpwZ6xUdTLOIEGZmZnWth/PNYUrX5hPZXUNRWVVvHXjsYSF2YVm07CoiDCuyOrJv77K5ofcQruN2fzIl1NJrwDvACcCiMhQEXnZh+VzAe9RYNKB7XXmyQSmiMgm4ELgKRE5z4dthKSyympufX0JOwrLyOgUx2vXjbFuLkyT3HBib5LaRfLPGevcjmKCiC+FIUxVPwWqAVR1BTDUh+UXAP1EJENEooCLganeM6hqhqr2UtVeeIrQzar6gQ/bCDlrdhZxxfPfszZvP/+5ZhRvTxrLqF7JbscyrURCTCTXjM3gy9X5rN5R5HYcEyR8KQzbRSQD5/SPc3dSk69qqmoVcCueu41WA2+p6koRmSQi1lDuCKgqd7+9nHV5xTxw7mDGDejsdiTTCl09thfx0RE8+XW221FMkIsXqmYAABbtSURBVGjyNQY8A/U8B6SKyDXABGCFLxtT1WnAtDqv1XuhWVWv9mXdoUJVqapR/ue9H2gXFc4P2wp55IJh/HxUD7ejmVYqsV0kl2f15N/fbuC3+0rtLjbj0xGD4ikGt+Fp+TwTuMIfoUzDHp+RTb/ffco7i3J5ee5m0pJi+ekx1nDNNM+FI9NQhRlr7C5x41theF9Vq1T1Hafl81PA0f4KZg5VUl7Fv7/dAMCwtETCw4TbTulLVIQv/4zGHKpPSjw9O7bjq9V5bkcxQaDRU0kichEwAkgQkUHAOlU92JH7ZGC4H/MZR2FpJe8vzuVARTXv3TyWET06sKu4nI5xUW5HM22AiDBhSCrPf7eRDQXF9EmxHvVDWVO+as4GVuEZse0xYL2ILBaRj4FSf4YzHuvy9jP24Rk88NEqjkpPZESPDgB0io+2TvFMi7n+xN7ERIbzzy/Xux3FuKzRIwZV3Qa8LCIbVHU2gIgk42kBvcbP+UJe4YFKrn95IeFhQlR4GLeM96nfQmOarFN8NOce1Y2pS7dRUVVjpyhDmC8D9cz2erwH2OOXROZHizbv4dbXl7CruJwpN2QxLC3J/rMavzp5YGfemL+FhZv2MLZvJ7fjGJcc0UA9eG5TXQ6sUNVyP2ULaYu37OWqFxaQkhDN69dnMbKnNVoz/je2T0eiI8J4a+FWKwwhzAbqCSKqyv6ySpbn7uOq5+fTKT6KKTdkWUtmEzBx0RH84vgMPli63UZ4C2G+NHDLVtX3ncdv+yNMKKuuUX7zznI+WLqNdpHhJMVF8vr1WXRpH+N2NBNibh3fl7cWbOWxL9bxyrVj3I5jXBCwgXpMwzbuKuGpr7N5d3Eu/TrH0zE+itevy7IWqMYVcdERTDqpD7PW72L+RruUGIp8OWIYgqfTvHtEZBGwFFiqqnb00Awl5VWc8/gsSiqq6d0pjmm3nYAIdhuqcdXlWT2ZPCuHx75Yy5QbjnU7jgmwRo8YROQpEbkJeBRPt9jeA/XYcWYzFB6o5PnvNlJSUU1Gpzj+9NNhhIWJFQXjutiocG4e14d5OXuYYyO8hZymHDEsxdO6+WI8Rw3FwA947krydbAeg+ci8/SVeUx6dREAZw5N5enLR7qcypjaLhndg3/PzOGfM9bbHUohpikN3CZ7PxeRdDyFYhhwNp6xm00TrdhWyG1vLCFnVwkAA1MT+ONPhricyphDxUSG84vje/HnaWtYsa2QoWk2wluo8OUaAwCqmotnNLZpjc1r/ktVeWnOJv48bQ0d46O487T+XDAy3S4wm6D288we/N+X63lx9iYevegot+OYAPG5MBjfrMvbz9dr8vlmbQFzc3ZzysDO/P1nR9HBOr8zrUBiu0guGJHOmwu2cufp/e2LTIiw/hX8RFV55LM1nP6Pb3n40zVs3XuA+88dzHNXZVpRMK3KL47PQATO/dd37CwsczuOCQA7Ymhhn63YwacrdpK7t5RFm/dy8aju3Hn6AFISot2OZswRyegUx3s3j2XiE7N5+pts/jjRl6HeTWtkhaEFfbRsO798YwkJMRFEhAm3ndyXX53W324/Na3ekG6JXDgynTfmb+WmcX1JTbQW+W2ZFYZmqqquYW7Obj5bsZO3Fm5lZM8OvHbdGGIiw92OZkyLumV8X95ZlMszMzfwgN1J16ZZYTgCW3Yf4Ks1eZzYP4X7p65k1vpdxEWFMyajIw+fP8yKgmmTuie344IR6bw+fws3jetj/Xi1YVYYfLTvQAUTn/yOvQcqf3zt92cP4vKsnlYQTJt3y/i+vLM4l2e/zeH35wx2O47xEysMPnhzwRbuefcHAJ6+bARLtu7j6O5JnDWsq8vJjAmMHh3bMWFoKm8vyuWuMwbYl6E2ygrDYZRVVjNl/hZ2FJaRntyOP32yik7xUdx+Sj/OHNaVM60gmBB0RVZPPlm+g9veWMITl46wUQXbICsM9SivqmbK/K08+XU2+fvLiQgTqmqUvp3jee26MXZu1YS0MRnJ3H3GAP42fS0fL9/O+SPS3Y5kWpgVBi97Sir4Zm0+//oqm427ShjdK5nHLzmGY3okkVNQQp+UePt2ZEKeiHDzuD68v2Qbz87ayNnDuxIdYaeU2pKQLQzZ+ft5Y/5W8orKGJ2RzPSVO1meW8j+sioSYyN58ZpRjOuf8mMbhEFd27uc2JjgISLccWo/bn19CXe+tYwnLh3hdiTTgkK2MGzZc4BX5m0mKTaSj5fvoEv7aE4b1IXzR6QzqGsCHeOtpbIxh3PO8G5sLCjh0S/WceHIfMYN6Ox2JNNCRFXdztAsmZmZunCh78NCVFTVUFVTQ2xkONn5xXRNiiU+OmTrpDFHpLyqmjP/OYuKqhqm33EicfZ/qNUQkUWqmlnftJA9YR4VEUa7qAhEhH5dEqwoGHMEoiPC+esFw9m2r5S/TV/rdhzTQgJaGERkgoisFZFsEbm3numXichy52eOiFgH8MYEucxeyVx1bC/+M2cTy7buczuOaQEBKwwiEg48CZwJDAYuEZG6TSc3Aiep6nDgIWAyxpigd9cZA+gYF8WfPllNeVW123FMMwXyiGE0kK2qOapagWdI0IneM6jqHFXd6zydB9gN0sa0AvHREdxz5kDmb9rDpc9+T16RjdvQmgWyMKQBW72e5zqvNeRa4NP6JojIDSKyUEQWFhQUtGBEY8yRuiizO09cegwrtxcy7m/fsGjzHrcjmSMUyMJQ36AE9d4SJSLj8RSGe+qbrqqTVTVTVTNTUlJaMKIxpjnOGd6Nz24/kZjIMP4zZ7PbccwRCmRhyAW6ez1PB7bXnUlEhgPPARNVdXeAshljWkivTnGcNawrX6zayfZ9pW7HMUcgkIVhAdBPRDJEJAq4GJjqPYOI9ADeA65Q1XUBzGaMaUGXjO6Bqmec6KKyysYXMEElYIVBVauAW4HpwGrgLVVdKSKTRGSSM9sfgI7AUyKyVER8b7lmjHHd0LRE3rghi90lFTwwdSXF5VVuRzI+CNmWz8YY/7v5tUVM+2En/bvE88b1WdbVTBCxls/GGFc8fvExTL5iJJt2H+DXby2jpqZ1fxENFVYYjDF+ExEexulDUrnvnMHMXFfAI9PXuB3JNIF1EGSM8bvLx/Rg3c79/HtmDl3bx3D1cRluRzKHYUcMxhi/ExEe+MkQTh3UhT9NW82q7UVuRzKHYYXBGBMQ4WHCXy8cTmJsJOc9NZsbX1nIruJyt2OZelhhMMYETHJcFB/eejyXju7B12sL+NWbS63TvSBkhcEYE1BpSbE88JMhPDRxCLPW7+KmVxfb3UpBxgqDMcYVPx/Vg/vPHcxXa/K58Jk5ZOcXux3JOKwwGGNcc/XYXlx7fAaLt+zjDx+ucDuOcVhhMMa4RkS475zB/P7sQczZsJvj/vIV//PeD+wotM733GTtGIwxrrt6bC+iI8KYl7OH9xbnMjt7Fy/9YjQZneLcjhaS7IjBGOO6iPAwrji2F09eNoI3bzyW/P1lnPzoN3y4dJvb0UKSFQZjTFA5unsSX981jlE9k7njzaWMeOgL7np7GZXVNW5HCxlWGIwxQadrYizPXpXJHaf0p29KPO8symXYA9N5d1Gu29FCgnW7bYwJep+v3Mlz321k/sY9nNQ/hd+dPYj+XRLcjtWqHa7bbSsMxphWobyqmpfnbOaJr7M5UFFF+5hIzjsmjbSkWDJ7dWB4epLbEVuVwxUGuyvJGNMqREeEc/2JvTl/RBqPz1jP1r2lPP/dRgAiw4W7Th/AdSf0JjxMXE7a+llhMMa0Kh3jo/njxKEAlJRXsau4nIc+XsXDn64hd28pPTu2o1N8NGcOSyU6ItzltK2TnUoyxrQJv35rKe8t/u/trZ3io/nrhcM4eWAXF1MFL7vGYIxp88oqq5m7YTeDurZnff5+Hp62hlU7isjoFEdKfDSjMjpwRVYv2sdGEBsZjkhon3KywmCMCTkHKqp4cfYmVmwrZHthGSu2FVLt9OKa3iGW845OY2yfjlSrktW7I5HhoXX3vhUGY0zIyykoZuqy7UQ5XW98t76Ag719d06I5rTBXbj+hN6kJEQDEBfdti/BWmEwxpg68ovKWJ5bSFWN8t7iXL7L3kVpZTWqntHmLsrszm2n9GVdXjErthUyMDWBoWmJdHYKR2s/FWWFwRhjGrF9XylvL8wlPAzyisp5Zd7mQ+YRgZ7J7QgT4aqxveieHEvB/nK6JcWSlhRLWodYVKGyuoaEmEgX3kXTWWEwxhgfrc/bz8x1BVRWKxdlprNxVwnvLMplzobdhAls2n3gkGUyOsVRUl4FwF2nD2DR5r0AnDG0CyN7JLN17wF6dYojPghOU1lhMMaYFqSq5O8vJ3dvKSnx0XyXvYvC0kq+XpNPQXE5G3eVANA+JgJV2F9eRXiYUF2jxESGMX5AZxZt3ktsVDiDUttzoLKamIgwbjypN31S4kmIifR7Qz0rDMYYE0A7CkvZvPsAR3dPokaVpVv28dnKnfTtHM/izXtZvGUfA1MTiAwPY+nWfdSoUlpZzb4DlQB0jIsiNTGGwtJK+nb2FIqTB6bwxao8IsLC6NmxHQkxEQxPTyKrd8cjymiFwRhjglxxeRWfr9zJ3gOVfJ+zmwMV1STERJCdX0xeURlFZVV0jIsiNiqc7ftKqVG48cTe/M9Zg45oe9ZXkjHGBLn46AjOH5EOwLXHZ9SaVlJeRU5BCX06x9EuKoKKqhoqqmuIjfRPlx9WGIwxJsjFRUcwLD3xx+dREWFERfivQV5oNfUzxhjTqIAWBhGZICJrRSRbRO6tZ7qIyOPO9OUiMiKQ+YwxxgSwMIhIOPAkcCYwGLhERAbXme1MoJ/zcwPwdKDyGWOM8QjkEcNoIFtVc1S1ApgCTKwzz0TgZfWYBySJSNcAZjTGmJAXyMKQBmz1ep7rvObrPIjIDSKyUEQWFhQUtHhQY4wJZYEsDPU146vbiKIp86Cqk1U1U1UzU1JSWiScMcYYj0AWhlygu9fzdGD7EcxjjDHGjwJZGBYA/UQkQ0SigIuBqXXmmQpc6dydlAUUquqOAGY0xpiQF7AGbqpaJSK3AtOBcOAFVV0pIpOc6c8A04CzgGzgAHBNoPIZY4zxCGjLZ1WdhufD3/u1Z7weK3BLIDMZY4ypzVo+G2OMqcUKgzHGmFqsMBhjjKnFCoMxxphaWv1APSJSABw6andw6ATscjvEYVi+5rF8zRPs+SD4MzYnX09VrbeFcKsvDMFMRBY2NEJSMLB8zWP5mifY80HwZ/RXPjuVZIwxphYrDMYYY2qxwuBfk90O0AjL1zyWr3mCPR8Ef0a/5LNrDMYYY2qxIwZjjDG1WGEwxhhTixWGZhKRF0QkX0RWNDB9nIgUishS5+cPQZZPRORxEckWkeUiMiKQ+ZwME0RkrZPh3nqmu70PG8vn+j50ciSLyBcist753aGB+TaJyA/OvlwYwHw/E5GVIlIjIg3eYtnY/g6SjG7tw7+JyBrn7+x9EUlqYL7m7UNVtZ9m/AAnAiOAFQ1MHwd8HMT5zgI+xTN6XhbwfYDzhQMbgN5AFLAMGBws+7CJ+Vzdh145/grc6zy+F3ikgfk2AZ1cyDcIGAB8A2Qe6f52O6PL+/B0IMJ5/Eh9/8YtsQ/tiKGZVPVbYI/bORrShHwTgZfVYx6QJCJdA5MOgNFAtqrmqGoFMMXJFCyaks/tfeid4yXn8UvAeS5kaJCqrlbVtY3M5urfQxMzukZVP1fVKufpPDyjXNbV7H1ohSEwjhWRZSLyqYgMcTtMHWnAVq/nuc5rwbZ9t/ZhU/K5vQ8P6qLOiIfO784NzKfA5yKySERuCFi6pgmWfdmYYNiHv8BzpFpXs/dhQAfqCVGL8fRJUiwiZwEfAP1czuRN6nktkPcwN2X7bu7DpuQL2D4UkS+B1Hom/c6H1RynqttFpDPwhYiscY4s/ZpPVT9syirqea1F92ULZASX96GI/A6oAl6rbxX1vObTPrTC4GeqWuT1eJqIPCUinVQ1WDrmygW6ez1PB7YH0/Zd3odN2T8B24eqempD00QkT0S6quoO51RWfgPr2O78zheR9/GcemiRD7XD5Wsiv+/LFsjo6j4UkauAc4BT1LmoUEez96GdSvIzEUkVEXEej8azz3e7m6qWqcCVzp01WUDhwdMRAbIA6CciGSISBVzsZPqRy/uw0Xy4vw+9c1zlPL4KOOTbr4jEiUjCwcd4LmbWe8eaS5qyv13l5j4UkQnAPcBPVPVAA7M1fx8G+qp6W/sB3gB2AJV4KvW1wCRgkjP9VmAlnjsD5gFjgyyfAE/iuYvhBw5zJ4YfM54FrHMy/M55LZj2YWP5XN+HTo6OwAxgvfM72Xm9GzDNedzb2Y/LnH36uwDm+6nzN1gO5AHT6+ZraH8HU0aX92E2nusHS52fZ/yxD61LDGOMMbXYqSRjjDG1WGEwxhhTixUGY4wxtVhhMMYYU4sVBmOMMbVYYTDGGFOLFQZjjDG1WGEwQUtEbhSRHU6f99ki8oGIRInIWBH5o4/rel5EznYe+7x8E/Id/BnW3PW2QK5eIlIqIksbma+viPxQ57VoEdkoIkOc91MhIp38m9gEG2vgZoKWiDwJLFTVF0UkDE9LzotV1eeBUURkETBRVXNbON9iVX2+pdbprDdcVaubsXwvPONXDG1kvjA8XbInq2qN89ovgZGqerXzfBOeltzB0reXCQA7YjDBbBiwxHncF0/XE+tE5G0ROR7AGcXqf0VklojsFJFTndf7i8h34hll61dA6sGi0MTlB4nIt+IZKetuEcluIF+938oPs94MEflQRBaKyHwRGeCV6TER+Rr4n4a2LyLDRGS213ZGiMhXje3I+rbrFIMtQC9nnljgTuCBxtZn2jYrDCaYDQFeFpHVwCLgGvX0tDoUT59EOI/3qeoJwM3AZSISAbwK/FpVh+HponuN13qbsvxrwO2qOhxP3zj1dZI2BHjR6zSSd7/89a03EnjOyZWJ5wP44LCLw4BiVR0P/OUw218J9BGRcOf5o8Bdh9uJjWx3NTDQeXwLMFVVNx1ufabts263TVASke5AvvPBiIhcCdwnIucCkapaKCLtgETgH85iEcA+4HxgtarOd15fCZQ664lp4vLLVPXg0coq6nRhXTdfnWkNrfc8PMXkXfF0FhsBzHIyJQMPOvM3uH1VrRGRlcAQEekHbFHVxY3sznq360xbDQwQkW/xFIasRtZlQoAVBhOshuP5QDxoGZ7THEO8Xh8CLPI6Hz8czzfr4XiOMA4aiWcMX3xY3vsU0VDgs3ryraF+Da33KDw9Xda6JiEiI/GME13lNf/htj8POA7PkciEBjJ4q3e7jtXAycDtwGuqmteE9Zk2zk4lmWA1DM+HFuL5mnsV8KXz+nJnnqHU/gAd7kzb7Uw7+KF7idd8TV2+v7P80cDleApT3XwNFYaG1rsDOMO56HvweoHUyUQTtj8P+F/gfVXd1kAGbw1tFzz7eDSeYSL/1oR1mRBgRwwmWA0DTnJuMVU8H4Z3AQ8B33vN873XMkPxfDPfCExzbtdci+c0zup6ljnc8p+IyAJgLrBJVXMayHem81yBE1S1+DDrXQGMB1aLSCmwQlUvd25xne81/yuNbH8NnvECHjlkr9Xvhfq260xb6+T9naoWNnF9po2z21WNqUNE4p0PeETkbiBRVX8fLNsXkSeABar6UgPL96IJt6s2Mcsm7HbVkGOnkow51K9EZKVzxNELz1GK69sXkT4isgaIbagoOKqBRGmkgdvhiEiss3wkUHOk6zGtkx0xGGOMqcWOGIwxxtRihcEYY0wtVhiMMcbUYoXBGGNMLVYYjDHG1GKFwRhjTC1WGIwxxtTy/5s2ytXJKaEsAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEaCAYAAAAVJPDdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeXhU5fXA8e+Z7HvIJOySoKKAu+KOK7SKa7XVVgMitkUSa2mrv9Y2bW1t6b7RVkDqUoFpXdpq0bpb6woq7iKKiGyyJ5B9m8n5/XFvxplJAhlIJjPJ+TxPnszc+87MyU0yZ95dVBVjjDGmnaevAzDGGBNfLDEYY4wJY4nBGGNMGEsMxhhjwlhiMMYYE8YSgzHGmDCWGEzcE8ddIrJLRF4RkdNE5IOQ8+tEZHIM4jhVRD4UkToR+Vxvv15fE5EzRWTTfjz+eyJye0/GFAv783Pv7zWLF5YYYsB942oUkVoR2S0iL4nILBHxhJT5q4i0uG867V83h9yuFxGNOD+qm69/hYi8F3HsyS6O3SQi3xWRRyLOfdjFsS9Ff0WiNhH4DDBSVU9Q1edV9dAYvG6kW4A/q2q2qj7o/j4O7oM44k5nb4iq+jNV/UpfxWT2nSWG2LlQVXOAYuAXwHeAOyLK/Mp902n/+nH7beAwt0x+yPkNoQ8WkatF5K+dvPazwDgRKXLLJQNHAZkRx04GnnO/ThWRJPfcUCAFODbi2MFu2f3ivvaeFAPrVLV+f19rPxUDK3viibrxMxvTZywxxJiqVqvqUuCLwHQROTwGr7kZWAuc7h46FucN7tmIYx5gBfAqTiI42j13OvAM8EHEsY/c5w4jIiXup+mZIrJZRLaIyA0h538kIv8QkSUiUgNcLSLDRWSpiFSJyBoR+apb9svA7cDJbi3px3uqrouIx631fCQilSJyn4gUdOc6icgJIrLMrdVtEZE/i0iqe+4j4EDgITeOZe7D3nLvf9Etd4GIvBlSMzwy5PnXich3RORtoL6z5CAih7k1tyoR2SYi33OPp4nIH9zrudm9neaeO1NENonIDSKy3Y19hnvuJBHZ2p7Q3WOXuDHs8Xk7iS2shuTWcn8qIlnAo8DwkNrscPf3vCSk/EUistK9Nv8TkXER1+ZGEXlbRKpF5F4RSe/m7+08EXlPnBr5JyJyY8i5i93fR437N3Gue3yGiKxyH7NWRK7dw/MPF5F/isgOEflYRL4eci7DvQ67xKmBH9+dmOOdJYY+oqqvAJuA02L0ks/xaRI4HXgeeCHi2HJVbVHVFuDlbpTfW23hLGAM8FngJgnvB7gY+AeQD/iAv+Ncj+HAF4CficgkVb0DmAUsc2tJN+/lNb8OfA44w32uXcCte3lMuwDwTaAQp/Y0CSgHUNWDgA04Nb9sVT3ZfcxR7v17ReRY4E7gWsAL3AYsjXijvQI4H6fm5w99cRHJAZ4CHnNjPxh42j1dAZyEk5iPAk4Avh/y8KFAHjAC+DJwq4gMUtXlQD1wdkjZK4G/dfN598qtyU0BNofUZsM+MIjIITi/428ARcAjOEk2NaTY5cC5wGjgSODqboZwB3CtWyM/HPiv+5onAIuA/8P5OzsdWOc+ZjtwAZALzAB+7/7+wojT3PsQ8BbOtZ0EfENEznGL3Awc5H6dA0zvZsxxzRJD39oMhH6avdH9NLVbRHb28GuF1g5Ow3mjfz7i2LP7Ub4zP1bVelV9B7gL502x3TJVfVBV23DeiCcC31HVJlV9E6eWMC26HxFw3pQrVHWTqjYDPwK+0Nmn80iq+pqqLldVv6quw3ljPyOK1/4qcJuqvqyqAVW9G2jGeeNt90dV3aiqjZ08/gJgq6r+1r0Otar6snuuFLhFVber6g7gx4Rfn1b3fKuqPgLUAe39MH/HvfZu8jnPPdad5+0pXwT+o6pPqmor8BsgAzglpMwfVXWzqlbhvBkf3cnzdKYVGC8iuaq6S1Vfd49/GbjTfc02Vf1EVd8HUNX/qOpH6ngWeILOP6QdDxSp6i3uh6a1wF+A9r61y4E5qlqlqhuBP3b3gsQzSwx9awRQFXL/N6qa734VducJRGReezIB5gFXhiSXt0OKPgccKSKDcN6olrn/JMPcYxMJrwE8B0x0zxWp6ofAS8Ap7rHD2XuNYWPI7fU4n4I7OzccqFLV2ojyI/by/J0pBh4IuSarcGoCQ/b2QBE5REQedpteaoCf4SStaF77hpDrvxs4gK5/7kgHAB91cW44zjVpF3k9KyNqIA1Atnv7b8Clbs3lUuB1VW1/rr09b08Jex33A8FGwn/HW0Nuh8a/N5/HSXbrReRZEWmvzXV5PUVkiogsd5vsdruP7+x3XYzTRBb6O/0en/49Dafj33nCs8TQR0TkeJx/ihf253lUtbw9meA0e/wtJLkcGVJuLU4NZSawQVXr3FPL3GPZwPKQp16G0zQxE3jRfY6akOfYrKof7yW8A0Juj3IfGwwp5PZmoMD9NBta/pO9PH9nNgJTQq5Bvqqmq2p3nms+8D4wRlVzcd4AJMrXnhPx2pmq+veQMntazngjTpNEZzbjvEm1i7yeXVLV93DesKYQ3owU7fM2AJkh94eGvsxewgh7HRERnL+Pffkdh1HVV1X1YmAw8CBwn3uq0+vpJsh/4tRahrj/O4/Q+e96I/BxxO80R1XPc89voePfecKzxBBjIpIrIhcA9wBL3GaWWHke+Jb7vd0L7rEVoc0b7u0VeyjfndFIPxCRTBE5DKcd997OCrlV8JeAn4tIujgdtl/G6XuI1gJgjogUA4hIkYhc3H7S7eS8uovH5gA1QJ2IjAXK9vJa23A6pNv9BZglIieKI0tEzo9IeHvyMDBURL7hdgrniMiJ7rm/A993f55C4IfAki6fqaO/4fS/nA7cH3I8mud9E6dGmuR24oY2s20DvCKS18Vj7wPOF5FJIpIC3IDTzPZSd4IXp+P7zE6Op4pIqYjkuU1UNTg1RHD6Hma4r+kRkRHu7zUVSAN2AH4RmYLTD9aZV4AacQYNZLg/++HuB7v2n+u7IjJIREYC13fn54l3lhhi5yERqcX5BFIB/A7nzTKWnsX5VBVaS3nePdbZG3205Tt7/BqcDtTfqOoTeyh7BVCC88nyAeBmVX2yG68RaS6wFHjCvd7LgRPBeRPB6RRe3sVjb8T5RF2L8ybfaSIL8SPgbreJ4XJVXYHTz/BnnE7vNXS/AxW3Ke0zwIU4zSof4nTgA/wUJ1G/DbwDvO4e666/A2cC/1XV0P6raJ53thvbbpy+iQdDYn/ffY217vUIa45S1Q+AqcCfgJ3u81zoDnTYI/cNt86NrzPTgHVu898s93XaB3jMAH4PVOP8PRa71/nrOG/qu3B+50s7e2JVDbixHg187MZ+O05tGpw+mfXuuSeAxXv7eRKBqG3UY3qYiJTg/KOkRI686UsiMhG4TlWv2GthEzdEZCpwmKp+t69jGSgsMZgeF6+JwRjTPdaUZIwxJozVGIwxxoSxGoMxxpgwCb+QV2FhoZaUlPR1GMYYk1Bee+21napa1Nm5hE8MJSUlrFixoq/DMMaYhCIiXc7StqYkY4wxYSwxGGOMCWOJwRhjTBhLDMYYY8JYYjDGGBMmZolBRO4UZ9vBd7s4LyLyR3G2dXy7s92UjOktPp+PkpISPB4PhYWFFBYW4vF4KCkpwefbl0VejUlcsawx/BVn276uTMHZBnIMznr/82MQkxmAQpNASUkJkydPZurUqaxfvx5VpbKyksrKSlSV9evXM3XqVLKzsy1BmAEjZvMYVPU5d3G1rlwMLFJnjY7lIpIvIsNUdUtMAjQDgs/nY+bMmTQ0NACwfv161q9fT3L+MNJGjCUpaxAkJROo20XLtjW0bnf2Iqqvr2f6dGc739LS0j6L35hYiKcJbiMI3yJvk3usQ2IQkZk4tQpGjeoXGyaZGKmoqAgmBUnNJPuoz5Jz9HmkFHS+m6W/die1bzxC7ev/IdBcz+zZsy0xmH4vnhJDZ9vqdbrCn6ouBBYCTJgwwVYBNN3i8/lYv349IGQf9VnyT7+KpMw8mja+S82rD9C8aSX+mh1owE9yjpe0kYeRNe40Bp1+FbnHf46qJ2+jctWzlJSUsGHDBkaNGsWcOXMsUZh+J54SwybC904dSTf3tDVmb8rLy5k/fz6ezHwKL/gWGaOPpWnju2y//2Zatq7pUN6/eyv+3Vupf/dpUgYfiPeccoou+j/qDjyO9Y/9Edz+h5kzZwLWvGT6l3garroUuModnXQSUG39C6YntCeF1CEHMezquaSNPIzKx/7Mtr/d1GlSiNS6fS1bl3yb3c8vIfvwsxnyxZ8iaVkANDQ0UFFR0ds/gjExFbMag4i07zlbKCKbgJuBFABVXQA8ApyHs09uA7HfD9n0Mz6fj1mzZlFXV0f6gcdRdPF3aWusYeuib9G6s8v1wzqnbVS/dA+tuzZTeP43GXzZj9h+7w/Q1ia3ecqY/iNmNQZVvUJVh6lqiqqOVNU7VHWBmxRQx3WqepCqHuFurG7MPikvL2fq1KnU1dWRceAEBl/yfVqrNrFl8Z6TQkpKCkuWLEFV8Xq9Hc43rHqOHUt/RdqwQyi6tAI8SQA2nNX0K/HUlGRMj2hvOgJILz6KoksqaNmxju33VNBWv7vLx3m9Xu66665gf8HcuXMR6TgmonH1Miof+xMZJccw6KxrAGc464wZMyw5mH7BEoPpV0KTQkphMUWXfI/WXZ+w/d7v09Zc36F8ampqsIawc+fOsE7k0tJSZs2a1enr1L/zFDWvPkjuhIvJGn8mAK2trdbfYPoFSwym3whNCknZBQy+7GbaWprYfv+POk0KxcXF3HnnnXscUTRv3jyWLFlCdnZ2h3O7/ncXTRvfpeCz5STnDQGw/gbTL1hiMP1CaFIgKYWiS3+AJy2b7f/4EYHanWFly8rKUFXWrVvXrWGmpaWl1NbWdkwQbQF2PvxbALznfxPEg4hYc5JJeJYYTMLz+XwsWLAgeL9g0kzSho1h58O/CS5p0c7r9TJv3rx9ep32BFFWVhY8FqjZQdWTC0g/4HByjr0AVWX27Nn79oMYEycsMZiEV1FRgbPEFmQddhY5x0yhevn9NK55JaxcZmYmc+fO3e/Xa29eale/8r80rn2N/NOmkpRdQGVlpdUaTEKzxGASXnu7for3AArOuY6mDe+w+7nFYWWKi4tZuHBhj81QLi0tpbi4OHi/6sn5SFIygyZ9FcA6oU1Cs8RgElrwk7knmcILb0Rbmti59FegbUD0/QnRmDNnTvC2f/dWqpfdR9bY00gbMZ7169dbrcEkLEsMJmH5fD6uucaZR5A/8UpShxxE5WN/IlC/K1hmX/sTuqO0tDRsnkPNqw/gr60k/yxn0r7NazCJyhKDSVizZ8+mpaWFtBHjyT3x89S+9TiNa14Ong9t6ukt7X0bANraTPWLfyN9xDgyxpxEa2urdUSbhGSJwSQcn89HTk4OlZWVSHIa3vO/gb9mB7v+e3tYudCmnt4SmXzq3n6S1sqNDDpjOojHOqJNQrLEYBKKz+fj6quvpq6uDoC8iVeSMmg4lY/MRVsag+W8Xm9MlsKeM2dO+LIZ2sauZ+8mxXsA2UdMBqwj2iQeSwwmoVx77bX4/X4AUoccRO7xn6P2rcdp3vhOsExKSkqPDEvtjs6WzWj8cDnNm1eTe9JlIB6bDW0SjiUGkzB8Ph/19e7SFuKh4NzrCTRUs/uZO4NlPB5P2EJ4sdA+r8Hj+fTfqXr5faQMGkbm2NMAZ2a2MYnCEoNJGKFNMrnHf460oQez68kFwXWQRIRFixb1yW5qpaWlLFq0KNis1Pjhy7TsXE/eyZcBwoIFC6yvwSQMSwwmYWzYsAGA5Lwh5E0spWH1MhpWvxQ8P2vWrD7dYrO0tDRklJJSs+x+UotKyDj4eFTV+hpMwrDEYBJGQUEBgDO7WNuoeurT9ZHKysp6dc5Cd4WOUqpf9Rz+3VvJO/lywFZeNYnDEoNJCOXl5VRWVpJ+4AQyx5xE9Yv3EKitBOInKUDEEFlto/qVf5E2fCxpI8aG9UEYE8/sL9XEveDqqUnJFEyaSWvlJmpW/BvYv9VSe0NkU1b9u0/T1lRHznEX0dbWZv0MJiFYYjBxzefzMX36dFSV3OMvIaVgOFVP3QZtzpDVqqqqPo6wo9DmJG1tpvbtJ8g89FSScrw2E9okBEsMJm6Vl5czbdo0AoEASTlF5J38Reo/eJGmdW8Ey4waNaoPI+xc5Izrutf/A0DO0efZTGiTECwxmLjUviNb+yifQWdfAyJhy16ISEyWvYhWaWkpXq83eN9fvY3GNa+QffS5SHKq1RpM3LPEYOJO5I5s6aOOJGvsadQsv59AzY7g8b4enronkTOva15bSlJmHpnjzrBag4l7Ero6ZCKaMGGCrlixoq/DMD2opKTk06Gd4mHY9D/gScvkk9vLINAaLBfvf7uFhYVUVlYG7w/78q1oSyNbF99IcXEx69at67vgzIAnIq+p6oTOzlmNwcSd0PH+2UdMInXIgex69q9hSSEWS2rvr8haQ91bT5A2fCwphcU2p8HENUsMJq74fL7gshKSmkHeadNo2vQeDe+/ECwTr30LkSL7GupXPoMGWsk+6hzA1k8y8csSg4krs2fPDjYR5Z34eZKzCzp0OMdz30KkuXPnBhNdW2MNDauXkXXYWZCUYusnmbhlicHEjfbZzQBJOUXkHH8J9Sv/R8uW1YDTfLR48eK4mtC2N+HrJ0HdW4+TlJFD5iEn2/pJJm4l93UAxkDHkUiDzpgOwK5n7wZI6M7a4uJP+xSa1r+Nf/dWso88h4ZVzwUXBjQmnliNwcSFioqK4Cfr1GGHkHXYmdS++gCBWmd4aiL0KXQlfJc3pe7tJ8koOYrk/KFxOUHPGEsMJi6EjtIZdPZX8NdVUb38H0DstunsLZG7vNW9+xTaFiDvmCkJnfBM/xXTxCAi54rIByKyRkRu6uR8nog8JCJvichKEZkRy/hM3wjtgM0cO5H0keOpfn4x2toEdBz2mYhOPfVUUlJSAAjUVtK49jXSx53JtbNmWQe0iTsxSwwikgTcCkwBxgNXiMj4iGLXAe+p6lHAmcBvRSQ1VjGa2PP5fFxzzTXOnaQUBp1xNS3b1lL3ztPBMolcW2hXUVFBa+un8zDq332a5BwvgcKDmTFjhiUHE1diWWM4AVijqmtVtQW4B7g4oowCOeI0yGYDVYA/hjGaGJs9ezYtLS0A5B53Ecn5Q9n1zB2gbUBiTGTrjshO5oY1r9DWVEfW+LNobW210UkmrsQyMYwANobc3+QeC/VnYBywGXgHmK3qvkOEEJGZIrJCRFbs2LEj8rRJED6fLzg81ZORS94pl9Ow5hWa1r8VLNNf2uA7dDIHWqn/4EUyDz0FSUmzmdAmrsQyMUgnxyIXuzkHeBMYDhwN/FlEcjs8SHWhqk5Q1QlFRUU9H6mJidBVRvMnXomkpLPrmTuDxxK90zlU+MgkR/3KZ/CkZpBx8EmIiDUnmbgRy8SwCTgg5P5InJpBqBnAv9SxBvgYGBuj+EwMhU5mS/aOJPvoKdS+8Sj+qk0ApKSk9ItO53aRI5MAmjeuxF+9nezDz7LJbiauxDIxvAqMEZHRbofyl4ClEWU2AJMARGQIcCiwNoYxmhjw+XzMnz8/eH/QmdegLY1Uv/g3ADweD3fddVe/qS20mzdvHkuWLAk5otSvfIb0kmPwZOZbc5KJGzFLDKrqB74GPA6sAu5T1ZUiMktE2j9K/QQ4RUTeAZ4GvqOqO2MVo4mN0Cak9OKjyDz4BKqX3UdbYw0AixYt6ndJoV1paWlYh3rde88gniSyxp9uzUkmbth+DCbmgm3t4mHY1XPxpGYE91rwer3s3Nm/Pwv4fD6mTZsWnOk99KrfA7B10TcTeukPk1hsPwYTN0I/EWcdPonUwaPD9lroT/0KXYlcWK9+5TOkDRtDsnekNSeZuGCJwcRUezOSpKSTf/o0mj5ZFdxrITs7u982IUUKbU6qf/85VNvIGnuaNSeZuGCJwcRM6Eik3BMv7bDXQujqqv1d6PDVtvrdNG9cSebYiTY6ycQFSwwmJsrLy4MjkZJyvOSecCn1q56jZfMHQP+as9AdHZqT3n+B1MJiUgpHWXOS6XOWGEyvi9xrIf+0aYh42PW/vwaPDYS+hUhJSUnB2w2rX0S1jcxDJ4YdN6YvWGIwvS50u87UIQeRfcRkalYsJVCzHRh4tYV2gUAgeLutfjfNG94lc+ypYceN6QuWGEyvCl0PCWDQWV8m0FBN9bL7AGfo6kCsLUDHBQLrP3Cak1ILi60D2vQpSwymV1177bXB2xkHn0B68ZHsfsGHtjQAMGvWrAFZW4CO6yc1rH4JbQuQceipYZMAjYk1Swym10yePJn6+nrnjieJQWdeQ2vlRurefAyAsrIy5s2b14cR9q3IDujQ0UmVlZVWazB9xhKD6RXl5eU8/fSnm+3kHD2FFO9Idv3X2WvB6/UO6KTQrkNz0vvPk1o4ipTCYqs1mD5jicH0uMhRSJ60LPImXknjujdpXOssXzJQ+xUiRe430bB6GdoWIHPsqVZrMH3GEoPpcRUVFWFNJLmnfBFPerZTW8BZPXWg9itEKi0txev1Bu+3NeymeeO7ZI09DcBqDaZPWGIwPS50G8vk/KHkHnsh9e88ReuOj4HwDmnTsfZU//4LpHgPIKWw2GoNpk9YYjA9xufzUVJSElZbyD/jarQtwO7nnX0IJk2aZH0LESJrDZ82J00EsCUyTMxZYjA9wufzMXPmzLDlHNJGHUHW2InUvPwPAnVVlJWV8dRTT/VhlPErtNbQ1rCbpg3vkuUmBlsiw8SaJQbTIyoqKmhoaPj0gHgomHwt/t1bqXnlXzYKaS861Bo++LQ5yVZcNbG218QgIgXd+MqPRbAmfoX2KwDkHHs+qUUlVP33DtTfQlVVVR9Fljjmzp0bnPDWPtmtfcVV64Q2sdSdGsNmYAXw2h6+3u6tAE388/l8YTN4PZl55E8spfHj12n8cBkAo0aN6qvwEkbohLe2huqw5iTrhDax1J3EsEpVD1TV0V19AZV7fRbTL7X3LbS1tQWP5Z9+FZKSTtVTCwFnPaTI8fqmc6ET3href95pTioqAawT2sROdxLDyT1UxvRDs2fPDutbSB06huwjP0PNa0vxV20CBvZ6SNEKTaANH4aPTrJOaBMre00MqtrUE2VM/xO5cioIBZ+5lrb6aqpf/DuAdTpHqbS09NOd3RqqadrwDlmHOonB9mkwsWKjksw+8fl8XHXVVWHHsg4/m7ThY9n17F1oSyNgS1/si9B5IA3vv0CKdyQpRSW2T4OJmf1KDCIytKcCMYnD5/Mxffr0sH4FT1oWg868muZP3qf+3WeAgbsBz/4K62dwRydljZ1ow1ZNzOxvjeGOHonCJJTZs2d3+PSaf8bVeDJyqXpyPqBkZmZabWEfhe7T0NZYQ9OGt8k81Bm2ah3QJhb2KzGo6vk9FYhJHOH9CpA6fCw5x0yh9rWHaNn2EV6vl4ULF1ptYR9F7tPwaXPSaOuANjFhfQwmKh2aMjxJeM+5Dn/NDna/4MPr9bJz505LCvspvDlpWbA5CZy9LozpTd1ODCLyqojcISLfEJGzRaSoNwMz8SlyBm7uhItIHTyaqqcWoi2N1nzUQzptTnITw4IFC6yvwfSqaGoMFwP3A6nALGCdiFi9dgApLy8Pa0ZKyi0i79RSGj58mcYPl1FWVmY1hR7SaXNSwQhSBo+2vgbT67qdGFR1s6o+pqq/UtXLgQnA7b0XmoknkbuyARRMngVA1VMLbL5CL+i0OelQm+xmel80TUlhi92o6irgsB6PyMSlyF3ZMsacTOaYE6l+0UegZoc1IfWCDs1J698i093ZzYaumt4UTVPSvSKySUSeF5F5IvI7YGw0LyYi54rIByKyRkRu6qLMmSLypoisFJFno3l+0zt8Pl/YJ1RJzaTgM9fSsv1jalYstfkKvaS0tJRZs2YF7zd88CIpBcOtOcn0umiakk5W1ZHADOBJYCVwQXcfLyJJwK3AFGA8cIWIjI8okw/MAy5S1cOAy7r7/KZ3tE9mCzXorBkkZQ2i8rE/IdpmtYVeFNo89+noJKfWELnUuTE9Jerhqqq6RlUfUNU7VHVTFA89AVijqmtVtQW4B6dDO9SVwL9UdYP7Wtujjc/0rMjJbOnFR5Fz9BRqXn2Qli2rbYG8GGjva/i0OcnpZygoKOjLsEw/Fst5DCOAjSH3N7nHQh0CDBKR/4nIayJyFZ0QkZkiskJEVuzYsaOXwjUQPplNUtLxnns9rVWfUP2C075tHc69b86cOaSmpgLu6KRBw0kZfCC7d++2fgbTK2K5VpJ0ckwj7icDxwHnA+cAPxCRQzo8SHWhqk5Q1QlFRTadojf4fD5ycnLCjuWffhVJeYOpfPSPqL8lbNSM6T2lpaWkpKQA4c1JgUDAdnYzvSKWayVtAg4IuT8SZ3e4yDKPqWq9qu4EngOO2r8QTbR8Ph/XXHMNdXV1wWNpIw8jd8JF1L72MM2bVgLY5jsxVF9fD0BbU21Yc1Lk8iTG9IRYrpX0KjBGREaLSCrwJWBpRJl/A6eJSLKIZAInAqv2J0YTvYqKClpaWoL3JTkN75Sv49+9ld3PLQKwyWx9qOH950kZNIzUIQf1dSimn4pmHoNXRMpEZIaInCAiGdG8kKr6ga8Bj+O82d+nqitFZJaIzHLLrAIew9lD+hXgdlV9N5rXMfsvcvJU/hnTSSkYQeVjf0Jbm2wyWx/wer3B2w2rl6EBP5njTgc6Wb/KmP0UTY3hAaAI+Bnwa6BaRN6P5sVU9RFVPURVD1LVOe6xBaq6IKTMr1V1vKoerqp/iOb5zf6LXKAtveQYcidcRM2KpTStfwsRseGpfSD0mrc11dG07s3gonrWz2B6WjSJIUdVbwG2qeoZwBXAXb0Tlom19s7m+fPnB4950rPxnvcNWnZuYPezfwVs/+a+EnnN699/juS8IaQOO4TKykqrNZgeFU1iaN/XuVlEMlT1n8B5vRCTibHOOpsBCj5bTlJmHpUP/xb1t7BkyRJrQupDYWsnffgy6m8ly21OslnQpidFkxh+IyIFwP1QpGYAACAASURBVL3AnSJyPR3nIZgENHv27LDOZoDMcWeQNe50dr/go2XbRxQXF1tNoY+FjgLT5noaP37dHZ0ktqie6VHRLInxT1WtUtXfAY/gDD2NnLlsEozP5+sw5DEppxDvZ8to2vQeNS//E7ChqfGgtLQ0vBP6/edJzikkbcRYW1TP9Kh9Gq6qqotV9duqurKnAzKx4/P5uOqqiMnl4qHwov8D8VD5n9+BtpGdnW21hTgxd+7c4IqrDWteRv0tZI49DVW1TmjTY/aaGETk9Z4oY+JLe79CW1tb2PH8iVeSPvIwKh+/Ff/urSQlJXXYh8H0ndANfLSlkcaPVjjNSeKxTmjTY7pTYxgnIm/v4esdoLC3AzU9K3ISGzgL5OWefDm1bz1Ow6pn8Xq93H333VZbiDOhndD17z9PcnYBaSOdhYqt1mB6QnI3ynRnz4XA3ouYeBLZWenJzKfwghtprdzIrqcWsmTJEksIcWrOnDlMnToVgMaPXqGttYmssafRvPHdYK3Bfndmf+y1xqCq67vxFc3y26aPRU5iA6HwwhuQtEx2/vtXFORZn0I8C+2E1tZmGj96lcxDTwVx/p1t6KrZX7FcdtvEgc72bs47+XIySo5h11O30bpzvc1sTgChv6OGVc+TlJVP+qgjANsP2uw/SwwDSPsopNC9m9NHH0veaaXUrXyGurefsG06E0RpaWlwdFLj2hW0NTcE94NOSkrqy9BMPxDNInpPiYgtgZ2gysvLmTp1atgopOS8IRRe+H+0bl9H1WN/tnWQEkxwdJK/hcY1r5B56CngSQ7bcc+YfRFNjeHbwO9F5C4RGdZbAZme5/P5wtZAAmcp7aJLvgci7HhgDupvtnWQEkzY6KRVz5KUkUvG6GPDJsEZsy+imfn8uqqeDTwMPCYiN0e79LaJPZ/Px7Rp0zocLzjnOlIGj2bnQ7/GX72NsrIyWwcpwYRu+dn48esEGqrJOuxMampqbD6D2S9R9TGI06j5ATAfuB74UEQ6vuuYuNDefBTapwCQc+wFZB9+NtUv/I2mta/Z/goJqrS09NPtV9sC1K96noyDT8QvyTafweyXaPoYXgA+AX6Ps3je1cCZwAkisrA3gjP7rrPRR+DsrzBo0ldp+HA51S/dS0pKivUrJLCqqqrg7fr3/ocnJY3MQ06xWdBmv0jkp8kuC4ocDqzUTh4gIqtUdVxPB9cdEyZM0BUrVvTFS8e1wsLCDovjpXgPYOi03+Cv3s5W37ehtYnFixdbv0ICKykpCRueOnzmX/BXb2X7vT/A6/Wyc+fOPozOxDMReU1VJ3R2LpqmpEWRSUFETnJvRrP3s+llna2Y6snIpegLN6OtLWz/xy1oS6MlhX4gctXb+veeIb34KJKyC6zWYPZZdxbRu1xEfgHkiMg4EQkdJL0QQFXX9laAJnrXXntt+IGkZIouqSA5u4Dt//oJgdodlJWVWVLoByKX4q5f+T9EPGSOOwOwtZPMvulOjeFF4D1gEPA7nA7n10XkYaCxN4Mz0SsvL6e+vj7smPec60k/4DB2PvIHUmo+sZ3Y+pnQPiL/rs00b15N9mFnAlitweyTaPoYTlXVF93bBcBo4H1Vrd/zI3uX9TGE83g8YaOQ8k+/iryTL2f380tI/uBJa3Pup0L7lHKOu5CCydey+Y5yWnduoLi4mHXr1vVtgCbu7Fcfg4jc7N5UEckBcHdye62vk4IJV15eHpYUco69gLyTL6f2zUepfukeG33Uj4X+butXPY+2BcgafyYAGzZs6KOoTKLqTlPSE+73bwAvi8hqEXlQRG4RkS/0Ymymm3w+Hzk5OWGzmzPHTmTQ5Jk0rF5G1RPz8Xg81qfQj4X2NbQ17KZp3RtuYhAKCgr6NDaTeLqz7PYy9/vlqjoeOAK4BfgIOGlPjzW9r30ntrq6uuCx9FFHUnj+DTRvWsXOh34N2taxQ9r0O3Pnzg3OhK5f+T+S8waTdsDh7N692/oZTFSimeD2rIjkqmozcAJOZ/T3ei0ys1ftM5tDd2JLHXowRZd+n9Zdm9nxz1tQfwvZ2dnW2TwAlJaWkpKSAkDD6mW0NdeTfcRkAoEAs2bN6uPoTCKJZh5DvqrWiMhxwFdxEsNfeicsszeTJ0/usDBeStFoBl/+EwKNNWy//4e0NdcjIrZn8wDSPiJN/c3Uv/ccmWNPRVIzqKurs1qD6bZoEkOriCQDVwG/VNWbgcN6JyyzJz6fj6effjrsWIr3AIZ88SdoaxPb/v49ArXOCBVbMXXgqnvnSTwp6WSNOx2wnd1M90WTGP4EvAVcADzkHsvu8YhMl3w+H4WFhcH9ftsl5w9j8Bd/imob2+6pIFCzHcBWTB2AQie7tWxZTcuO9WQf8RnARieZ7otm2e27gROBw1W1UUQOBpb1WmQmTHsnc+RSF8l5QxhyxRwkKZnt91Tg37UZj8djk9gGqMghyXVvP0HaiLGkeA+w0Umm26LpfD4bmAv8VERmAHmA9WjFyLXXXhvWyQyQPGg4Q678BZKSzrZ7v09r5UY8Hg+LFi2y5qMBqrS0lLKysuD9+pXPoAE/WUdMtn0aTLdF05S0BGeTnuXAgcAPgZW9EZQJN3ny5A7LXKR4D3CSQlIK2/7+PVq3f0xqaqolBcO8efM+ndPQWEPDmpfJPvxsWgNttnaS6ZZoEsMaVX1AVe9X1R+o6sWqenA0LyYi54rIByKyRkRu2kO540UkYBPouuhoLhrNkCt/Aaps/ftNtO5Yh9frpbm52ZKCASL2aXj7SZKyBpFx0PG2dpLplmgSw7Mi8k13F7eouauy3gpMAcYDV4jI+C7K/RJ4fF9epz/x+XxcddVVYcdShx3CkCt+hvqb2fa3m/BXbrLNdkwHo0aNCt5u/Ph1/LWV5Bx1LmArrpq9iyYxHAaUAVtE5D8iMkdELovi8Sfg1DrWqmoLcA9wcSflrgf+CWyP4rn7lfYlLqZOnUpbW1vweProYxnypTm0Ndez1XcT/t1b8Hq93HXXXVZTMGHC9mnQNureeoz0A48lOW8IlZWVlJeX911wJu5FMyrpUlU9BGdV1ZuBD4luSYwRwMaQ+5vcY0EiMgK4BNjjjCwRmSkiK0RkxY4dO6IIIf75fD6mT58etsQFQNb4Mxn8+R/i37WZrUtuJFCznbKyMnbu3GlJwXQQuU9D3VtPgCrZR08BYMGCBdakZLoUzaikMSJyB/BrVV2hqn9V1RuieK3OmqAi1/z+A/AdVQ3s6YlUdaGqTlDVCUVFRVGEEN98Ph/Tpk0jEAj/8XOOv4TCC2+kedNKtv7tJtrqdzNp0iQbjmr2KLR5MVBXScOHy8k+8jOQlIKq2oQ306VompIWA/8ATgdnD2gRWRTF4zcBB4TcHwlsjigzAbhHRNYBXwDmicjnoniNhNU+TyFsfwzxMOjsr1Bw9pepX/Uc2+6/GW1pxOv18tRTT/VdsCYhdKg1vPEISZl5ZB16KkDYXtHGhIomMXhU9VEgAKCq7wKHR/H4V4ExIjJaRFKBLwFLQwuo6mhVLVHVEpwkVK6qD0bxGgnH5/NRUlLSYTE8Sc2g6NLvk3v856hZ8W92Lv01BPzW0WyiEvq30rT+LVorN5J9rLNFu8cTzb+/GUii+cvYLCKjcZt/3NFJGd19sKr6ga/hjDZaBdynqitFZJaIDMiJcu2ro0Z+ckvKHczQqb8m48DjqHz8VnY9/RdAycrKso5mE5XIv5XaNx8lfcQ4UgaPpq2tzfoZTKei2dqzBLgdZ6hpBXAuTi0impFJPS5Rt/b0+Xwd1jwCSBsxlqJLvg9Jyex88Oc0rX8LsHWPzL4rKSkJfvjwpGUx4rq7qV/5DFWP34rX67XtXgeo/draM4TiJIOv48x8fhaYtv/hDUydjSXPOuxshnzp57Q1N7B18Q2WFEyPCB262tZcT/17z5J12Fl40rNtwpvpVDSJ4QFV9avqP9yZz/OAo3srsP6qfY5C2GJ4SckUfLacwgu+RfMn77F18Q34qz6xxfBMj4jshK5dsRRPSnpw6KpNeDOR9poYRORyEfkFkCMi49yZye0W9l5o/U97n0LoHIWk3CKGXvkrco45j+rl97Pt3h/Q1lSLiNi6R6bHhHZCt+5cT+Pa18g57kJISrZag+mgOzWGF4H3cHZs+x3woYi8LiIPA429GVx/Ul5e3mHHtfTRxzLs6rmkeEew/V8/Zfezd4M6M51tgx3TkyJrDTWvPkBydgFZ484AbBMfE26viUFVP1HVRcDFqjpFVQ8EJuPMfj67twPsDzokBfGQN/FKBl/2IwK1lWy5+5s0frgccDZaseYj0xvChq6ue5OW7R+Te7wzTcjmNJhQ0SyJ8WLI7SpVfU1V6/f0GNMxKSTnDWHIlb8g/9QrqX/3v2xdfCP+XZsREZYsWWJLXJhe01mtIXXwaNJLjgGw9ZNMUFQb9YjIHSLyWxGZISLHiUhabwaXqNo7mEUkLClkjT+TYTP+SGpRMTuW/orKR/6A+psBazoysTF37lzaF0iuX/Uc/roqck+4BID58+dbX4MBopvHsBm4DkgGjnS/Dot2T4aeFk/zGHw+H7Nnz+6w/aakZlDwmTKyDz+bpk0r2fnQb4P7MosIixcvtqRgYiZ05fzcky5j0BnT2Xzn9bTu+NjmNQwge5rHkBzF86xR1Qfc2/fvf1j9i8/nY+bMmTQ0NIQdTy85Bu+515OU42X380uoXnZfsIM5KSmJu+++25KCiani4uJgn0LdG4+Qd9IXyDvlcnb++5cdPtSYgSlmG/X0d7Nnzw5LCpKaScG51zPkiz9B/c1s9X2b6pfuCSaFrKwsSwqmT0ROeKt97WEyDz2VZO9IAGtOMjHdqKffmjx5ctgnrfQDJzD8K/PIPmIy1cvvZ/NdX6dl8wfB82VlZdTV1VlSMH2itLSUrKys4P2aFf9GW5vJO/mLgE14M92b4DZPRMqA3+Isix26Uc+JvRte/PL5fBQWFiIiwT2ZPZl5eM//FkMu+xFtTfVsXXyjMzch0Bp8nC1vYeLBbbfdFrzd1lhD7RuPkDXudJIHDbcd3ky3agxv4tQWfgasxVkZ9WZgDBAfvb4xFLrt5qe1BCH76CkM/+ptZI07jd0v3cOWu2fTsvXD4OOysrJsfoKJG50NXdWAn7yTnUYA2+FtYOvOBLeFqvo1VT1DVQuBicB8oAY4v7cDjBc+n4+MjIwOS1qkDD6QodN+g/ec62jd9hGb77qe6ueXQMAPQGZmJkuWLLGmIxN3Qie8tdXvpu7NR8k67GySvSNRVWtSGsCi3qlDVTep6iOq+ktVHRCrq7avcdTU1BQ85knPYdDkaxk2/fck5w1m50O/Yds9FfgrNwXLeL1eFi5caAnBxKXIWkP1svvQ1mbyT3P+ra1JaeCyLZy6UF5eTlJSUodJaniSyTn+cwy/9i/kHHMetW88yid/mUX9e/8Le/ykSZNsFrOJe6ET3toaa6h59UGyDj2V1KHO9CSb9DYwdXuCW7zqjQlukydPDnYoh8oYczKDzppByqDhNK5dwa5n7qR154YO5SZNmmR7MpuEEbpsi6RmMOLa22nZ9hHb7/shgE1666d6aqOefq2zUUbtUoeOYcgVP2fwpRVooJVt9/2Q7ff/qENSaF8Az5KCSSTz5s0LNilpSyPVy+4jY/SxpJc4263YpLeBZ8AmhtBEICIRo4wcKYXFFF1SwbDpvyelcBSVj9/Kljuvp+nj14NlUlNTWbJkCapqTUcmYYV2RNe+8R9ad21h0NlfAXHeIqw5aWAZkE1JPp+Pa665hpaWlk7PJxeMIP/UK8kcdxra0kjNKw84k4BawrefsDkJpj/Jzs6mvt5ZMDnjkJMZfEkFlY/fSt2bj1pzUj9kTUkhfD4fV111VadJITl/KN4psxn+5XlkHHwiNcv/wScLvkz1S/eEJYWysjJU1ZKC6VdCJ701rl5G0/q3yT9tKpKWZbu8DTADosbQ1aqn7VKKSsg76TIyx06EtgC1bzxC9fL7aWuoDivn9XqZO3euNReZfquwsDD4f5JSNJphM+ZSu2Ipu/57Ox6Px7ab7Ud6anXVhFNeXs6CBQvoKvmljRhH7kmXkXnwCbQ1N1Dz6oPUrvg3gbqqsHLWZGQGirlz5zJ16lQAWnd8TN2bj5Fz3IXUvftfWrevZcaMGQCWHPq5fldj8Pl8VFRUdL1VoXjIOGgCucdfQvqoIwg0VFP72kPUvvYQbc3hG9JZDcEMRKG1Bk9aFsO/sgB/zQ62LrkRtM36G/qJPdUYEj4xJCcnayAQ2Gs5Sc0g+4jPkHPchaQMGoa/Zgc1rzxA3duPo63NYWVtnwQzkPl8vmCtASBz3OkUXfRtqp5cQO3rDwOwZMkS+/9IcP06MYjIHn+A5Pxh5Bx3IdlHTMaTlknTpveofW0pDauXQVvHhGK1BGPCaw0Agy/7MWkjxrH5jjICtZVkZ2dTW1vbhxGa/TXwRiV5ksg85BQGX34LI9ylKxrWvMyWu7/JNt+3aXj/hbCk0D7KyOYiGOOYO3cuqampwftVT8wD8eCd4iysV1dXZ6OU+rF+VWNIzh9G9lGfJfvwySRlD8Jfs4O6t5+k7q3HOnQog9UOjNmTyCal7KOn4D3nOpvb0E/086Ykj2aOnUj2UeeQUXI02hagcc0r1L71uDND2d1KE5xZynfeeaclAmO6qUOT0uW3kDZiPFvuuh7/7i02Yi+B9eumpJTBoym6+DukDBrG7ucW88n8a9jxwBya1q4IJgURoaysjObmZksKxkQhdKkMgMpH50Kbn8ILbwBPEvPnz7elufuhmCYGETlXRD4QkTUiclMn50tF5G336yUROWqvT9rSwLZ7KvhkwVeoXnYvgbrwSWzFxcUsXrzYPtUYsw9KS0spKysL3g/UVlL52J9IGz6W/DOmA87S3NnZ2dbn0I/ELDGISBJwKzAFGA9cISLjI4p9DJyhqkcCPwEW7u15jzxoBI3r3mTJksUUFxcjIhQXFwcXtlu3bp3VEozZD6GrrwI0fPAita8/TN4Jl5Jx8AkA1NfXM336dEsO/UTM+hhE5GTgR6p6jnv/uwCq+vMuyg8C3lXVEXt63t7Yj8EYE87n8zFt2rRPVxFISmHo1F+TnDeErXd/A3/1NsD2bkgk8dLHMALYGHJ/k3usK18GHu3shIjMFJEVIrJix44dPRiiMaYzpaWlzJo169MDgVZ2/vsXIELR53+ApGYAtndDfxHLxCCdHOu0uiIiZ+Ekhu90dl5VF6rqBFWdUFRU1IMhGmO6Mm/evLD+Bv/urez89y9I8R5A4QU30P4vbp3RiS+WiWETcEDI/ZHA5shCInIkcDtwsaraxw9j4khkcmha9ya7nv4LmWNOCuuMtuSQ2GKZGF4FxojIaBFJBb4ELA0tICKjgH8B01R1dQxjM8Z0U2RyqH39YWrfeIS8k75AznEXATZSKdHFLDGoqh/4GvA4sAq4T1VXisgsEWlvvPwh4AXmicibImK9ysbEociRSlVPLqDhg5comDyTzHFnADZSKZEl/MxnG5VkTN/obKTSkMtvIW3EOHY88DMaP3oFwBbci1PxMirJGNOPdDZSafs/f0LLto8ouuS7ZBx8IuAsuGd9DonFEoMxZp9F9jdoSwPb7v2Bkxw+910yxpwMWId0orHEYIzZL50nhx/SsnUNRZ+7iewjPwtYh3QiscRgjNlvnSeH79P08Rt4p3ydvFO+BFiHdKKwxGCM6REdkkNrE9v/9RPq3nmK/NOm4r3gBiQ5jUAgEN43YeKOJQZjTI+JTA60Bah85A/sem4RWePPYOjUX5GUO5i6ujprVopjlhiMMT2qQ3IAapbdx/Z/3EJy3hCGTf896aOPpb6+nhkzZlhyiEOWGIwxPa6z5NC0dgVb7v4mgfpdDLn8FgZNmolfhalTp1rtIc5YYjDG9IrOkoN/9xa2LvoWNSuWkjvhIoZO/wOpQw6ivr6eqVOn2pDWOGGJwRjTazpLDupvYdfTC9l27/fxpGUx9KrfkX/WNUhqhs13iBOWGIwxvaqz5ADOyqxb7ryOurefJO+ESxn+lflkHnIK8+fPR0QsQfQhSwzGmF43b948lixZgscT/pbT1lRH1eN/ZsviG2lrqKHoku8x5Iqfkzr8UObPn09aWpr1PfQBSwzGmJgoLS1l0aJFJCcndzjXsvl9ttz9DSqfmEdKwUiGTfstRZdUoLlDmTp1KiJCYWGhJYkYsdVVjTEx5fP5mDVrFnV1dZ2el5R0cidcTO6Jn8eTlknD6mVUL7+fli3hW7R4vV7mzp1LaWlpLMLud/a0uqolBmNMnykvL2f+/PmdnvOk55Bz3IXkHHchSRk5NK1/m9rXH6ZhzcvQFgiWmzRpEk899VSsQu43LDEYY+LWnpIDODWI7KPPJfe4i0jOG4y/tpK6tx+n7q3HCdR+uvuv1SCiY4nBGBPXfD4fX/nKV2hqauq6kHjIOPA4co45n/QDjwVVmta/Rf17/6Nh9TK0pRGArKwsbrvtNksQe2GJwRiTEPbW/9AuOW8I2Ud+hszxZ5KSP5S21mYa17xMw+qXaFz7OtrSECxrNYnOWWIwxiSU7iYIgNThY8kafyZZYyeSlJWPBlpp2vgujR++TMOaVwjUbA8rLyKoKsXFxcyZM2fAJgxLDMaYhOTz+Zg9ezaVlZV7Lywe0oYfSsbBJ5I55kRSvAcA0Fq5iab1b9G04W2aNrxDW2NN2MNSU1O58847B1yCsMRgjOk3srOzqa+v32u55EHDyTjoeNKLjyL9gMPxpGUC0LJtLU2b3qNl8/s0b/4A/+4tHR6bmppKS0tL8H5/bI6yxGCM6Td8Ph8zZ86koaFh74XbeZJIHXow6aOOJL34KNKGHRJMFIGGapo3f0DL5g9o3vYRrds/JlDXjRoKiZ0wLDEYY/oVn89HRUUFGzZsoKCggOrqavx+f/efQDykeA8gbcRY0oYfSurwsaQWjgqeDjRU07J9Ha07PqZl+1patq/Dv+sTtLU56ljjtU/DEoMxpt8rLy/ntttuo62tbZ8eL6mZpA4eTerg0aQMHk1q0WhSiorxpKQFy/hrdtBatYnWyk/wV22iteoTWqs2EajZCezfe2msax+WGIwxA040I5u6JB6SC4aTWlhMSsFIkgtGkFIwkhTvCDxpWcFiGmjFX70Df802/NXb8VdvI+B+91dvJ1BXxf4mjkj7O1/DEoMxZsBrb35av359jzyfJyvfSRIFI0jOG0py3mCS84aQlDeY5OyCsLIaaMVfW0mgrpJAXRWBul0ht6vw11USqK0Km3/RrRg8HhYtWrRPycESgzHGRIhqKGyUJDmVpNwikvOGkJw7mOT8ISRle0nK9pKcU0BSdkFYjaNdW2sTgfrdtDVUE2isCfve1lhDoKGaQENN8LY211NcXMy6deuij9ESgzHG7FlkosjKct64uzM0dl9ISjpJ2QVhX8nZXjyZeSRl5jnfM3LxZObiSUnv9Dk04KetsYZDDxxFbnoyuRkp5KankJuR7H7v7L5TrignvcvE0HFhdGOMGYBKS0v32iQT2hzl8Xj2uaMbQFub8O/ajH/X5r2WleQ0PJm5TsLICP+eN3gEh5w8nppGP1X1LazbWU9Nk5+axlb8bfv2wd9qDMYYs59Ch8+OGjWK8847j/vuu69XmqlC7amPQVVpbA1Q0+inpqmVmsZW97uf6sZWrj51tDUlGWNMPNrXYbb9ZlSSiJwLzAWSgNtV9RcR58U9fx7QAFytqq/v6TktMRhjTPT2lBhitueziCQBtwJTgPHAFSIyPqLYFGCM+zUT6Hr3DmOMMb0iZokBOAFYo6prVbUFuAe4OKLMxcAidSwH8kVkWAxjNMaYAS+WiWEEsDHk/ib3WLRlEJGZIrJCRFbs2LGjxwM1xpiBLJaJQTo5FtnB0Z0yqOpCVZ2gqhOKiop6JDhjjDGOWCaGTcABIfdHApEDeLtTxhhjTC+KZWJ4FRgjIqNFJBX4ErA0osxS4CpxnARUq2rHXTSMMcb0mpjNfFZVv4h8DXgcZ7jqnaq6UkRmuecXAI/gDFVdgzNcdUas4jPGGOOI6ZIYqvoIzpt/6LEFIbcVuC6WMRljjAkXy6YkY4wxCcASgzHGmDCWGIwxxoSxxGCMMSZMwq+uKiI7gJ7Zq69rhcDOXn6NnmBx9iyLs2clQpyJECP0TJzFqtrpDOGETwyxICIrulqFMJ5YnD3L4uxZiRBnIsQIvR+nNSUZY4wJY4nBGGNMGEsM3bOwrwPoJouzZ1mcPSsR4kyEGKGX47Q+BmOMMWGsxmCMMSaMJQZjjDFhLDG4ROROEdkuIu92cf5MEakWkTfdrx/GOkY3jr3FKSLyRxFZIyJvi8ixsY4xJJZzReQDN5abOjnf59e0GzHGzfUMialARJ4UkQ/d74O6KLdORN5xr+2KGMd4mYisFJE2EelyWOXern9viyLOPruW7uv/WkTed/8GHxCR/C7K9cz1VFX7cvpZTgeOBd7t4vyZwMMJEOd5wKM4u+GdBLzcR3EmAR8BBwKpwFvA+Hi6pt2MMS6uZ0RMvwJucm/fBPyyi3LrgMI+inEccCjwP2DCvl7/eIizr6+l+/qfBZLd27/s7Hfek9fTagwuVX0OqOrrOPamG3FeDCxSx3IgX0SGxSa6MCcAa1R1raq2APe4scWT7sQYL9czMqa73dt3A5/rw1g6paqrVPWDvRTr87+RbsbZ51T1CVX1u3eX4+xuGanHrqclhuicLCJvicijInJYXwfThRHAxpD7m9xj8RpHX17T7sQYL9cz1BB1dzZ0vw/uopwCT4jIayIyM2bRdV88XtuuxNO1vAanFhupx65nTDfqSXCv46wtUici5wEPAmP6OKbOSCfH+mJMcnfi6Otr2p0Y++R6ishTwNBOTlVE8TSnqupmERkMPCki77s1zh6xpxhV9d/deYpOjvX4te2BOKGXryV0L04RqQD8J55lYgAABbpJREFUgK+zp+jk2D5dT0sM3aSqNSG3HxGReSJSqKrxtuDWJuCAkPsjgc3xGEccXNPuXKs+uZ6qOrmrcyKyTUSGqeoWt1lrexfPsdn9vl1EHsBpauixN7M9xdhNMbm2PRBnr19L97n3GKeITAcuACap26kQoceupzUldZOIDBURcW+fgHPtKvs2qk4tBa5yR9OcBFS3NzvE2KvAGBEZLSKpwJfc2ILi4JruNUbi53pGxjTdvT0d6PCpV0SyRCSn/TZO52WnI9n6UHeuf5+Lh2spIucC3wEuUtWGLor13PXsq172ePsC/g5sAVpxMu+XgVnALPf814CVOD39y4FT4jROAW7FGZ3wDnsYaRGDWM8DVruxVLjH4uqadiPGuLmeITF7gaeBD93vBe7x4cAj7u0D3ev6lnuNK2Ic4yXu32czsA14PDLGrq5/vMXZ19fSjWENTv/Bm+7Xgt68nrYkhjHGmDDWlGSMMSaMJQZjjDFhLDEYY4wJY4nBGGNMGEsMxhhjwlhiMMYYE8YSgzHGmDCWGEzcEpFrRWSLuwb+GhF5UERSReQUEflxlM91h4ic796O+vHdiK/964j9fd4eiKtERBpF5M29lDtYRN6JOJYmIh+LyGHuz9MiIoW9G7GJNzbBzcQtEbkVWKGqd4mIB2dG55dUNeqNUkTkNeBiVd3Uw/G9rqp39NRzus+bpKqB/Xh8Cc4+F4fvpZwHZwn3AlVtc49dDxynqle799fhzPaOtzXBTC+yGoOJZ0cAb7i3D8ZZnmK1iNwvIhMB3N2sfioiz4vIVhGZ7B4/REReEGfXrW8CQ9uTQjcfP05EnhNnx6z/E5E1XcTX6afyPTzvaBH5t4isEJFXROTQkJh+JyLPAN/t6vVF5P/bu38QKc4wjuPfH96BFwIH9oEDMQns7SoqIoiIaWKKgNiJoqSNhWi0UkHUIkHEJmVISKKlCAGDhVgo4j8MGDz3bPQQQlAQFAQb9bF434WZdWZ3tHG9+32q3Z15n+dlivfZ993Zd9qSrhTyrJR0cdiFrMqbi8FDYCqfMwH8ABweFs/mNxcGG2Ut4A9JXeAW8F2kHVmnSfsWkV8/jYj1wPfANkljwClgb0S0SVt5zxbiNml/GtgdER3SXjlVm6a1gN8Ky0jFffqr4o4Dv+R+rSYNwL3HL7aB5xGxEfhxQP4ZYKmkRfn9CWDfoIs4JG8X+DK/3gX8FRFzg+LZ/Odtt20kSfoMeJwHRiTtAA5J+hYYj4hnkj4BJoGTudkY8BTYAnQj4kb+fAZ4keMsbtj+dkT0Zit36dvaur9/fcfq4m4mFZMzSpvKjgGXc5+WAEfy+bX5I+K1pBmgJWkZ8DAi/hlyOSvz5mNd4AtJl0iFYe2QWLYAuDDYqOqQBsSe26Rljlbh8xZwq7Ae3yF9s+6QZhg9q0jP9OUd2heXiKaB8xX9m6VaXdzlpB0vS79JSFpFepb0y8L5g/JfA9aRZiKbavpQVJk36wJfAbuB0xHxqEE8m+e8lGSjqk0atFD6mrsTuJA//zefM015AO3kY0/ysd6gu7VwXtP2n+f2K4DtpMLU37+6wlAX93/g6/yjb+/3AvX1iQb5rwHHgLMR8V9NH4rq8kK6xmtIj4s83iCWLQCeMdioagMb8i2mQRoM9wFHgeuFc64X2kyTvpk/AP7Ot2veIy3jdCvaDGp/TtJN4CowFxH3a/r3TX4fwPqIeD4g7h1gI9CV9AK4ExHb8y2uNwrn/zkk/yzp+QE/vXXVqv1alTcfu5f7eyAinjWMZ/Ocb1c16yPp0zzAI2k/MBkRB0clv6SfgZsR8XtN+yka3K7asC9z+HbVBcdLSWZv2yNpJs84pkizlA+eX9JSSbPARF1RyF4BkxryB7dBJE3k9uPA6/eNYx8nzxjMzKzEMwYzMytxYTAzsxIXBjMzK3FhMDOzEhcGMzMrcWEwM7MSFwYzMyt5A8pU7JVZCmgJAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plotXPS(energyXPS, XPS)\n",
    "\n",
    "energyGrid, dos, occStates, unoccStates, dosScaled, occStatesScaled, unoccStatesScaled, energyGridJdos, jdos, convGrid, finalPeakPreconv, finalPeakNorm, convGridShifted, finalPeakScaled = dosToPeak(energyGrid,dos,paramsDict,fermiLevel=fermiLevel,visualizeInput=visualizeInput, visualizeOutputs=visualizeOutputs)\n",
    "\n",
    "plotXPSFit(energyXPS, XPS, convGridShifted, finalPeakScaled)\n",
    "\n",
    "writeFiles(energyGrid, dos, occStates, unoccStates, dosScaled, occStatesScaled, unoccStatesScaled, energyGridJdos, jdos, convGrid, finalPeakPreconv, finalPeakNorm, convGridShifted, finalPeakScaled, paramsDict, name = name)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.7"
  },
  "nteract": {
   "version": "0.22.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
